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Abstract 

We study the dominant late-time behaviors of massive scalar fields in static 
and spherically symmetric spacetimes. Considering the field evolution in the 
far zone where the gravitational field is weak, we show under which conditions 
the massive held oscillates with an amplitude that decays slowly as t -5 / 6 at 
very late times, as previously found in (say) the Schwarzschild case. Our 
conclusion is that this long-lived oscillating tail is generally observed at time¬ 
like infinity in black hole spacetimes, while it may not be able to survive if 
the central object is a normal star. We also discuss that such a remarkable 
backscattering effect is absent for the held near the null cone at larger spatial 
distances. 

PACS numbers: 04.20.Ex, 04.70.Bw 


Typeset using REVTJrjX 


* Email: hir oko@ allegro. phys. nagoya-u .ac.jp 
fEmail: atomi@allegro.phys.nagoya-u.ac.jp 

1 



I. INTRODUCTION 


One of the most remarkable features of wave dynamics in curved spacetimes is tails. 
Scalar, electromagnetic and gravitational fields in curved spacetimes do not, in general, 
propagate entirely along the null cone, but are accompanied by ’tails’ which propagate in 
the interior of the null cone. This implies that at late times waves do not cut off sharply 
but rather die off in tails. 

In particular, it has been well established that the late-time evolution of massless scalar 
fields propagating in black-hole spacetimes is dominated by an inverse power-law behavior, 
as was first analyzed by Price [|I|. In a brilliant work, Leaver []2j demonstrated that late¬ 
time tails can be associated with the existence of a branch cut in the Green’s function 
for the wave propagation problem. Gundlach et.al. || showed that power-law tails also 
characterize the late-time evolution of radiative fields at future null infinity, while the decay 
rate is different from that of timelike infinity. Furthermore, it has been shown that power- 
law tails are a genuine feature of gravitational collapse ||-Q]: Late-time tails develop even 
when no horizon is present in the background, which means that power-law tails should 
be present in perturbations of stars, or after the implosion and subsequent explosion of a 
massless held which does not result in black hole formation. The existence of these tails was 
demonstrated in full non-linear numerical simulations of the spherically symmetric collapse 
of a self-gravitational massless scalar held. Gundlach et ah fH obtained the power-law tails 
for a massless held in fully nonlinear simulations at hxed r, Marsa and Choptuik § found 
them both at hxed r and along the event horizon, and Burko and Ori [pj at hxed r, the 
event horizon, and future null infinity to very high numerical accuracy. 

When the scalar held has a non-zero mass, the tail behaviors are quite different from 
massless ones. For example, as is well known, the tails exist even in Minkowski spacetimes, 
which is related to the fact that different frequencies forming a massive wave packet have 
different phase velocities [^]. If the background spacetime is curved, it is expected that 
interesting features peculiar to massive helds develop through the scattering due to the 
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spacetime curvature. 

The important role of massive scalar fields has been revealed in elementary particle 
physics. For example, in higher-dimensional theories the Fourier modes of a massless scalar 
field behave like massive fields known as Kaluza-Klein modes, and the recent development 
of Kaluza-Klein idea (e.g., the Randall-Sundrum model || in the string theory), strongly 
motivates us to understand the evolutional features due to the field mass in detail. In 
addition, scalar fields are astrophysically important, if boson stars made up of self-gravitating 
scalar fields prove to be viable candidates for dark matter ||. If such astrophysical objects 
become unstable and collapse to form a black hole, both gravitational and scalar fields would 
be presumably radiated. Further, as one example of curious effects peculiar to massive 
fields, it has been argued that unstable quasi-normal modes can exist 00- Then the 
time evolution of massive scalar fields in curved spacetimes, (in particular, in black-hole 
spacetimes) would become an important problem to be solved. 

Recently it was pointed out that the late-time tails of massive scalar fields in Reissner- 
Nordstrom spacetime are quite different from massless fields in the existence of the interme¬ 
diate late time tails |T^](see also |TB|). If the Compton wave length m -1 of a massive field 
is much longer than the horizon radius of a black hole with the mass M, namely mM -C 1, 
each multiple moment 0 of the field evolves into the oscillatory inverse power-law behavior 


-0 ~ t 1 3/2 sin (mi), 


( 1 ) 


at intermediate late times. It is clear from (|I[) that massive fields decay slower than massless 
ones, and waves with peculiar frequency u quite close to m mainly contribute to massive tail, 
while the dominant contribution to massless tails should be evaluated in the zero-frequency 
limit u> —> 0. Though the oscillatory power-law form (jl]) has been numerically verified at 
intermediate late times, mM <C mt 1 /(mM) 2 , it should be noted that the intermediate 
tails are not the final asymptotic behaviors; Another wave pattern can dominate at very 
late times, when it still remains very difficult to determine numerically the exact decay rate 
T2l . Pl ■ In the previous paper Jl4|, we have analytically found that the transition from the 
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intermediate behavior to the asymptotic one occurs in nearly extreme Reissner-Nordstrom 
background. The oscillatory inverse power-law behavior of the dominant asymptotic tail is 
approximately given by 


ip ~ t 5//6 sin (mt), 


( 2 ) 


independently of the multiple moment /, and the decay becomes slower than the intermediate 
ones. Then, the similar result for the decay rate has been obtained by considering massive 
scalar fields in Schwarzschild background (in the limited cases that mM C 1 or mM 3> 1, 
where M is the black-hole mass) [|15|] and massive Dirac fields in Kerr-Newman backgrounds 
T|. Asymptotic behaviors of massive scalar fields in dilaton black-hole backgrounds have 


been also discussed 17 


These results given in |fL4j [TT)|| suggest that massive fields in black hole backgrounds decay 
as f -5 / 6 generally at very late times. So it is an interesting subject to study how universally 
such a slowly decaying tail develops. It has been numerically shown Q that a power-law 
tail develops even when the collapsing massless scalar held fails to produce a black hole. 
This is an evidence for the late-time tail to be a direct consequence of wave scattering in 
far distant regions. In this paper we prove that the decay law f -5 / 6 of massive scalar fields 
can be essentially determined by the analysis in the far zone where the gravitational held 
is weak. However, we can also derive the conditions for the tails with the decay rate of 
f~ 5 / 6 to dominate as an asymptotic behavior. Considering the physical interpretation of 
the conditions, we can claim that any spherically symmetric black holes generate the same 
asymptotic tails, while the conditions may not be satishecl if the central object is a normal 
star. 

In Sec. [II] we introduce the Green’s function analysis to investigate the time evolution 
of a massive scalar held in any static, spherically symmetric spacetimes. In Sec. Ill we 
consider the approximation valid in the far zone, and we find the conditions for the tail with 
the decay rate of t -5 / 6 to develop. The final section is devoted to discussion, which contains 
a comment that the tail behavior of t -5 / 6 breaks down as the region becomes close to the 
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light cone. We discuss that the tail with the decay rate of t 5 / 6 can develop also in rotating 
black hole spacetimes. 


II. GREEN’S FUNCTION ANALYSIS 


A. Massive scalar fields in spherically symmetric spacetimes 


We consider the evolution of a massive scalar field in a static spherically symmetric 
background with the asymptotically flat metric given by 

ds 2 = — f(r)dt 2 + h(r)dr 2 + r 2 (d0 2 + sin 2 Odtp 2 ). (3) 

Here we do not assume the metric to be a solution of the vacuum or clectrovac Einstein 
equations. The scalar field <f> with the mass m satisfies the wave equation 

□<f> = m 2< f>. (4) 


Resolving the field into spherical harmonics 


$ 


r 




(5) 


hereafter we omit the index l of %) l for simplicity, and we obtain a wave equation for each 
multiple moment 


d t 2 


drl 


+ V{r) 




0, 


( 6 ) 


where r* is the Wheeler tortoise coordinate defined by 


dr * 
dr 



and V is the effective potential 


V = f 



1(1 + 1 ) , 2 
+ m 


(7) 


( 8 ) 
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The time evolution of the radial function 0 is given by 


= / [G(r*, r'; t)A{K, 0) + G*(r*, r'; t)0(r', 0)]dr' 


( 9 ) 


for t > 0, where the retarded Green’s function G is defined as 

G{n, r*; t) = 5(t)S(n ~ r*). 


d 2 3 2 t ; 

-h U 

dt 2 dr 2 


( 10 ) 


The causality condition requires that G(r*, r'; f) = 0 for t < 0. In order to obtain G(r*, r'; £), 
we use the Fourier transform 

G(r*,r';w) = J G(r*, r'; t)e luJt dt, (11) 

which is analytic in the upper half u; plane. The corresponding inversion formula is 

-i poo+ic 

G(r *, r*; t) = —— / G(r*, r'; (12) 

where c is some positive constant. Now the Fourier component of the Green’s function 
G(r*, r'; u) is expressed in terms of two linearly independent solutions for the homogeneous 
equation 


B 2 

W 2+uj2 ~ v 


A = o 


i = 1 , 2 . 


(13) 


The boundary condition for the basic solution ijji is that it should be well behaved on the 
event horizon if the central object is a black hole, and at r = 0 otherwise. On the other 
hand, the other basic basic solution 0 2 is required to be well behaved at spatial infinity, 
r —> oo. Using these two solutions, G(r*, r'; u;) can be written by 


G(r*,r*;u;) = 


A(A,u)A{r*,u) , r* > r' 

W ( u ) i ^(r*,w)$ 2 «,w) , n<r(, 


(14) 


where W(u>) is the Wronskian dehned by 


W(cj) = </>l<fe, r , - <Pl, r ,tp2- 


(15) 


The integrand in (O) has branch points at u — ±m. Considering the branch points, 
one may change the integration path in (|l^). First, if r* — > t, the path is closed in the 
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upper half of the u plane for the integration to converge. Since the integrand would have 
no singularities in the upper half plane, we obtain G(r*, r'; t) — 0 according to the causality 
postulate. If r* — r' < t, on the other hand, the path can be deformed to the curve shown 
in Fig. [[]. As will be shown later, the late-time tails are generated owing to the existence of 
a branch cut (in ip 2 ) placed along the interval — m < uj <m. 


B. The analysis in a region far from the gravitational source 


It has been found in the previous papers p~4|-|T7|j that the oscillatory power-law tails 
of massive scalar fields whose decay rate is t ~ 5 / 6 dominate at asymptotically late times in 
black-hole spacetimes. In this paper we show that the decay law can be simply derived by 
considering wave modes only in a far distant region, as a generic behaviors observed in any 
black hole spacetime. 

For that purpose, we assume 


r 

— » 1 , 
M 


(16) 


where M is the gravitational mass of a background field. Then, the expansion of the metric 
functions / and h as a power series in M/r leads to 


/=1 _^ + <£ + 0 ( r -) 




( 17 ) 


and 


A = 1 + ^ + e2 +0(r -3), 




(18) 


where M', Q and Q' are some parameters characterizing the background field in more details 
in addition to the gravitational mass M. Expanding (p~3|) in the same manner and neglecting 
terms of order 0[(M/r) 3 } and higher, we obtain the approximated form 

d 2 ^ 


dr 2 


U ^ = 0, 


(19) 


where 
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U = (m 2 - a; 2 ) 


2 2 Mu 2 2 M'(m 2 - w 2 ) A 2 - \ 


+ 


( 20 ) 


The coefficient A in (p0|) depends on the multiple moment l and the other parameters M, 
M ', Q and Q'. For example, in the case of Reissner-Nordstrom background with mass M 
and charge Q, we have 


A = 


l + — ) + 4m 2 M 2 — 12 u 2 M 2 — m 2 Q 2 + 2u 2 Q 2 . 


( 21 ) 


We keep the term of order of 0(M 2 /r ~ 2 ) in (|19|), in order to confirm that the decay rate 
of asymptotic timelike tails found in [|h|-[T7| is independent of A. Introducing the variable 
defined as 


x = 2 wr, 


where 


w = 


V- 




— a r 


(H) is rewritten by 


where k is 


d 2 1 k 
dx 2 4 x 


A 2 - 1/4' 


or 


■0 = 0 , 


Mm 2 

k = - 

w 


(.M + M')w. 


( 22 ) 


(23) 


(24) 


(25) 


III. TIMELIKE ASYMPTOTIC TAIL OF MASSIVE SCALAR FIELDS 

A. The wave modes 

Our aim is now to show that the tail with the power-law decay of is a generic 
asymptotic behavior in black-hole spacetimes by using the wave modes satisfying ©■ One 
may claim that the inner boundary condition to determine ijj± is missed if the analysis is 
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limited to the range (pj6|) . Hence, we treat -0i as a general solution for (|19D and reveal a 
condition which allows the excitation of the asymptotic power-law tail. Fortunately we will 
be able to prove that such a condition is always satisfied if the event horizon exist in the 
background spacetime. 

First, let us give 0 2 , by requiring that it damps exponentially for |u;| < m and to be 
purely outgoing for |u;| > m at spatial infinity. The outer boundary condition leads to 


<fc(ro,r) = Wk,i(x), 


(26) 


where W K> \(x) is the Whittaker function |H|, and the branch of w is chosen to be 


a/ m 2 — uj 2 , 00 < m, 


w = 


■i\J'co 2 — m 2 , a i > m. 


(27) 


Note that W Kt \{x ) is a many-valued function of w, and there is a cut in i/j 2 . The late-time 
tail is generated by the contribution from the branch cut in -0 2 , while -0 1 is a one-valued 
function of xzr, as was shown in ||T4| , |T5|| . This is because the late-time tail is a consequence 
of backscattering. On the other hand, we give 0 1 , using the Whittaker functions M K ^\ and 
free parameters a and b as follows, 


t/>l = aM KtX (x) + bMx-\{x), 


(28) 


where a and b will be determined if the inner boundary condition for 0 X is specified. Nev¬ 
ertheless the relation 


i/(i(ro) = ^i(e i7r ro) 


(29) 


should be required, since 0i is a one-valued function for w. Then, we rewrite the parameters 
a and b as 


a 


= jw 


1/2-A 


and 


(30) 
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b = fcro~ 1/2+A , 


(31) 


where j and k are some one-valued functions for w. In the following we will clarify un¬ 
der what kind of conditions for a and b the power-law tail with the decay rate of t ~ 5 / 6 
asymptotically dominates. 


B. Branch cut integration at timelike asymptotic regions 


As was shown in flT2| ; p~4| , |T5| , late-time tails are derived by the integral of G(r*, r(; uj) 
around the branch cut in Fig. [Ij. Using (|12|), (^6[) and (|28|), the branch cut contribution to 
the Green’s function is given by 



K;t) 



02 (r,w) 
W(w) 



ap + + 
ap + — bp _ 


•02 (r, e m w) 
W[e m zu) 

aq + + 6g_ 
aq + — bq_ 


e~ iuJt du 


AjkX 


e~ iujt dcu, 


(32) 


where 


and 


r(±2A) 
r(i± a-k) 


(33) 


Q± 


r(±2A) 
r(i± a + k) 


e ii r(i^A)_ 


Note that at very late times 


(34) 


mt 1, (35) 

the rapidly oscillating term e~ lujt leads to a mutual cancellation between the positive and the 
negative parts of the integrand (|32D, except the case that the other terms of the integrand 
also change rapidly with u. In fact, it is easy to see that neither oscillates rapidly nor 
changes exponentially with uj in the region 


ur -C k. 


(36) 
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Then, if k remains small, the effective contribution to the integral in ([32]) is claimed to be 
limited to the range \u — m\ = 0(1/1) or equivalently w = 0(^m/t) (see |I^Jh|Jl5[| ), and 
the intermediate tails becomes dominant at late times in the range 

1 


when the integral 


rnM < rut < 

(mM y 

should be estimated under the condition 
m 2 M 


(37) 


K ~ 


\Jm 2 — 


tx) A 


= 0(mM\fmt ) -C 1. 


(38) 


As was discussed in |14|,|T5|1, the small value of n represents that the backscattering due 
to the spacetime curvature is not effective at intermediate late times. It is obvious that 
the intermediate tails given by ([I]) dominate at intermediate late times (|37j), which was 
numerically supported by [0,0]. 


As was also discussed in , however, the intermediate tails cannot be an asymptotic 

behavior, and the long-term evolution from the intermediate behavior to the final one should 
occur. The asymptotic tail becomes dominant at very late times such that 

1 


mt 


m 2 M 2 ’ 


(39) 


when the effective contribution to the integral 

m 2 M 


k ~ 


arises from the region 

> A 


(40) 


\JrnY — uj 2 

which means the backscattering effect due to the curvature-induced potential dominates. In 
the limit of k —> oo, the term ( ap + + bpJ)/(ap + — bp _) becomes 

ap + + bp _ p + e l ' KK + 7 + e“ 


ap + — bp _ p_e~ mK + y_e J7rK ’ 

which includes very rapid oscillations as e ±lirK , and we have 


(41) 


rj± = T(2\)aK~ x e~ i7rX ± T(-2X)bK A e lnA , 


XAnX 


(42) 


and 
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7 ± = r(2A)<m _A e i7rA ± T(-2\)bK X e~ inX . 


(43) 


Such rapidly oscillatory behaviors are not seen in the other term ( aq + + bq_)/(aq + — bqJ), 
which is given by 

aq + + ar(2\)n~ x e~ mX + bT(—2X)n x e mX . . 

aq + ~ bq~ or(2A)K _A e _ * 7rA — bT(—2X)n x e inX 

even in the limit of n —> oo. 

Now we revisit the procedure through which the branch cut integration (|3^) leads to the 
tail with the decay rate of f~ 5 / 6 , as was shown in [14, HJ. For example, in the case of small 
mass field (mM -C 1) in a Schwarzschild background with mass M, r]± and become 


V± = 


r(2A) 2 r(l - 4 iuM) 
r(l/2 + A - 2iujM) 2 


(4m 2 M 2 ) _A e _ ” rA 


± 


r(-2A) 2 r(l - AiuM) 
T(l/2 - A - 2 iwMf 


(4 m 2 M 2 ) x e™ x (45) 


and 


7± = 


r(2A) 2 r(l - ^M) (imVyXe „ x ± r(-2A)V(i-4wM) (4m2M2)V , TA 


r(l /2 + A - 2 iujM) 2 y ' r(l/2 - A - 2 iuM) 2 

in the limit of n —> oo (see (39) and (40) in |T5|). Then the following inequalities 


(46) 


toil % It±I 


(47) 


are satisfied when uj 0 respectively (see (43) in |15]). When the inequalities ( [471 ) are 
hold, the term ( ap + + bp~)/(ap + — 6p_) can be expressed as the product of the rapidly 
oscillation term e ±2mK by e^, where tp remains in the range 0 < <p < 2ir even if hi becomes 
very large. The integrand (^) includes rapidly oscillating terms of e 2MTK and e~ lut , which 
means physically that the scalar waves have multiple phases owing to the backscattering by 
the spacetime curvature, and the contribution from these waves are canceled by those with 
the inverse phase, unless the phase of oscillation becomes stationary, i.e., 

— (ut =F 2irn) = 0 (48) 

auj 

for u ^ 0 respectively. We denote u satisfying ([48| ) by uq. Then, particular waves with the 
frequency ujq remain without cancellation, and contribute dominantly to the tail behaviors 
(see Fig. ^j). In the limit of |u;| —> m we obtain the solutions of (|48|) as 
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t ~ 


for lo ^ 0 respectively, or equivalently 


2TTLV 0 m 2 M 

: (m 2 -a;2)3/2 


(49) 


/ 2nM \ 1/3 


= \/ m 2 - uj ~ m I —— J 


(50) 


Approximating the integration 
obtain 


by the contribution from the close vicinity of cuo, we 


G c (r*,r';t) ~ — (mM) 1 ^ 3 (mt) 5//6 sin (mi + 0)'0i(r*, m)7/q(r(, m). (51) 

4v3Aj£: 

Thus we can confirm that the decay law of t~ 5 / 6 is a result of wave evolution in far distant 
regions r M. We also find the existence of a phase shift (ft given by 


4 > — - <p(wq) + )tt. 


(52) 


which modulates the basic oscillation with the period of 2k/ m (see |I5[). The multiple mo¬ 
ment l and metric component h{r) can affect only this modulation term in the asymptotically 
late-time evolution. 

Note that if the conditions ([47]) break down, i.e., either 


M $ Till 


(53) 


for u ^ 0, respectively, or 


\v± \ = h 


±i 


(54) 


are satisfied, all the contribution from scalar waves will be canceled more effectively. Then 
the tail with the decay rate of cannot survive. Therefore, the conditions are 
necessary for the tail to dominate. The physical implication is given in the next subsection. 


C. Physical interpretation of the condition for the tail generation 

It is easy to see that the inequalities ( |T7| ) are satisfied for the background spacetimes 
discussed in the previous papers (see (63) in [|b|], and (43) and (74) in [|l^]), in which the 
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tail with the decay rate of f -5//6 can dominate at very late times. Now we give the physical 
interpretation of ([|7]) to discuss what background spacetime allows the development of the 
late-time tail. 


Because A 2 in (|20|) is real, A should be either real or purely imaginary. From the expres¬ 


sion (PI) in Reissner-Nordstrom background, we find that the small mass ( mM -C 1) gives 
a real A, while the large mass ( mM 1) gives an imaginary A. The two inequalities (f47j) 
for u ^ 0 can be reduced to 


-{ab* -a*b) ^ 0 
A 


(55) 


when A is real, and 


7 


: (| 6 | 2 '- |«| 2 ) § 0 


(56) 


when A(= * 7 ) is purely imaginary, respectively. 

One may claim that the condition ( f47| ) changes according to the value of A. However, 
we can give the unified interpretation, independent of the value of A, by paying attention to 
-01 in the region 

1 


c or 


K 


(57) 


in addition to (fiG|). When all of the conditions (|1(|), ( |36p and ( |57| ) are satisfied, the second 
term which represents the Newtonian part in the effective potential (^0|) becomes dominant, 
compared with the other terms. In this region ipi is approximated by 


7 ~ Ae‘ 2 ^ + Be~‘ 2 ^*, 


(58) 


where A and B are 


A = 7T- 1/2 K' 1/, x 1/4 (2\)e-‘’ r/4 {n.~ x a r(2A)e"‘ rt - K A f>r(-2A)e i *-'} 


(59) 


and 


B = 7r -1 / 2 K -1/4 ;r 1/4 (2A)e - ” r/4 {/W A ar(2A)e i7rA - ^bY{-2X)e~ inX } , 


(60) 
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respectively (see Fig. |3|). In this region, independently of A, the mode clearly shows a wave 
behavior with the amplitudes |Al| and \B\ corresponding to the outgoing and ingoing parts 
for u > 0, while |Al| and \B\ correspond to the ingoing and outgoing parts for u < 0. The 
difference between \A\ 2 and \B\ 2 is 

\B\ 2 — |Al| 2 = 2i\7i(ab* — a*b). (61) 

when A is real. The conditions (|5^) under which the rapid oscillation of e 247rK survives are 
equivalent with the inequalities 


|B| g |U| (62) 

for u ^ 0 respectively. On the other hand, when A(= iy) is purely imaginary, we have 

\B\ 2 - \A\ 2 = |2A| 2 |T(2A)| 2 (|6| 2 - |a| 2 )(e“ 2i7rA - e 2 “ A ) 

= 87r 7 (|6| 2 -|a| 2 ), (63) 

which means the conditions ([>(]) are also equivalent with the inequalities (^2j). Therefore, it 
is sufficient to consider the inequalities ( p2| ) independently of A, as the conditions for the tail 
with the decay rate of t" 5 / 6 to dominate at late times. (|5^) means the amplitude of ingoing 
waves for is larger than that of outgoing waves, in the region where (0), (|36|) and (|57|) 
are all satisfied. 

The origin of the slowly decaying tail as of a massive scalar held can be considered 
a resonance by cooperation between dispersion and backscattering. It is a common feature 
when the scalar held has a nonzero mass that in far distant regions the effective potential (|20|) 
is a monotonously increasing function with r and the radial mode shows a wave behavior. 
If the central object is a black hole, the conditions (]6^) are surely satished because of the 
existence of the event horizon. So, we can conclude that this long-lived oscillating tail is 
generally observed in arbitrary spherical symmetric black-hole spacetimes. 
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IV. DISCUSSION 


We have found that whether the tail with the decay rate of t -5 / 6 develops at very late 
times can be judged relevantly by wave modes only in far distant regions. Then, even when 
the central object is a rotating black hole, only the parameters M' and A in the effective 
potential ([?!]) will be changed in far distant regions. Strictly, background spacetimes in this 
paper are limited in the class of static and spherically symmetric. However, since these are 
not relevant to the conditions (pn^), the same tail behaviors are expected to dominate also 
in Kerr spacetimes. 


We compare our analytical result with their numerical simulation |12|. As far as the 


intermediate late-time behavior is concerned, our result agrees with |12[. However they 
claimed ”SI perturbation fields decay at late times slower than any power law” in ||12|| , 


which disagree with our present result and previous ones [|14],[n| that the late-time tail of a 
massive scalar field is a power law with index —5/6. We believe that the integration time 
in 


121 is too short to find the true asymptotic behavior. 


Now we remark that the region of spacetimes where the tail with the decay rate of f _5//6 
dominates is limited. This feature can be understood by considering the behavior of -0In 
the region 


roo r > k(uj 0 ), 


(64) 


0i is reduced to 


-01 - \l —A 


hi 


2 ew 


^zur—K, In r (^ Ai tk 


(r/_e 4 ™ + + 


K 


2ew 


g— tot+k In r ^ ^ — iirn^, 


(65) 


and the saddle point at u = (|50| ) is disappeared because of the terms of e ’ wr ~ Klnr which 

change exponentially. Therefore it is obvious that the U 5 / 6 tail dominates only within the 
region 


r < M 1/3 t 2/3 . 


( 66 ) 
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What kind of behaviors dominates in the region r M 1//3 t 2 / 3 , in particular, near the 


null cone r —> t? Now we must hnd saddle points as solutions of the following equation 

d 


7 — (— iut + wr — k lnr) = 0 , 

ou 


(67) 


instead of (f|^). In general, solutions of ([T?|) are complex functions of t and r. However we 
can find a simple asymptotic solution as 

imt 


Uq ~ 


Vt [ 


2 _ 772 ’ 


( 68 ) 


for high frequency uq ^ m, which is compatible with the limit of 


f —> t, 


(69) 


where r is 


f = r + (M + M') lnr, 


(70) 


which is modihed due to red shift. The expression (^) means that we can calculate Green’s 
function using the saddle point integration by deforming the integration contour into the 
straight line AOB in Fig. [I]. When uq is a large value, considering the region 2uq r' 1 
also, Green’s function is reduced to 

/ p-iut 

(e-” R + e” R )^, (71) 

and we obtain 


imt 


uq ~ 


Vt 2 - R 2 


as the saddle point rather than (|68D, where R is 


R = r — r . 


(72) 


(73) 


Approximating the integration (|T^) by the contribution from the immediate vicinity of uq, 
we obtain 
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G(r*,r';t) 


0 i{—ut-\-y/ u 2 —m 2 R) 


‘li'Juj 2 — m 2 

CJ=CJi 

= i (- Wl t+ <v /^5Z^ fl ) 3/41/2 


do; exp 


a 2 


9a; 2 


ujt + \/u> 2 — m 2 Iij 


(u-UJxf 


LO—UJ 1 


T' ill mT lr \t + fi)~ 1/4 (i - fl) _1/4 


(74) 


Radial part of the scalar field <f> near null cone R ~ t, together with the geometrical factor 
1/r, behaves as 


where u is 




(75) 


u = t — R. 


(76) 


This behavior (|75|) is similar to the case of Minkowski spacetimes, except for f including red 
shift factor, instead of r. Massive fields near null cone decay more rapidly than f -5 / 6 . 

Finally we comment about late-time tail behaviors when the central object is a normal 
star like as a neutron star or a boson star. If the expression of ipi in (|28|) is assumed to 
be extended to the region r < M, then we must require Vh to be regular at r = 0. This 
leads to the equality |R| = \B\ which means that the amplitude of outgoing is equivalent 
with that of ingoing. Then the tail with the decay rate of f~ 5 / 6 never develops. Though this 
extension of ( |2§| ) may not be valid, we can expect the equality |A| = \B\ to be valid, unless 
some absorption of waves occurs in the inner region. This is a future problem to be checked 
by giving a background gravitational field with a regular center. 
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FIGURES 



to the power-law tail. 



integrand includes rapidly oscillating terms of e 2l7rK and e lujt , but the phase of oscillation is 
stationary at ui = ujq. 
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in the far region where all of the conditions 


( 0 ), 


@) and (H) are satisfied. 
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